Does ohmic heating influence the flow field in thin-layer electrodeposition? 
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In thin-layer electrodeposition the dissipated electrical energy leads to a substantial heating of 
the ion solution. We measured the resulting temperature field by means of an infrared camera. The 
properties of the temperature field correspond closely with the development of the concentration 
field. In particular we find, that the thermal gradients at the electrodes act like a weak additional 
driving force to the convection rolls driven by concentration gradients. 

PACS numbers: 81.15.Pq, 87.63. Hg, 47.27.Tc 



I. INTRODUCTION 

The electrochemical deposition of metals from aque- 
ous solutions in quasi-two-dimensional geometries has 
proven to be a valuable test bed to examine concepts 
of interfacial growth like fractal growth jl], g], morpho- 
logical transitions ||, Q, |5j or dendritic growth jg, 
The properties of the evolving deposit are in many cases 
sensitive to the presence of convection currents in the 
solution H, §, |§ ||, H, |||. There are two well- 
known origins of convection rolls: a) on small scales, 
spatial inhomogeneities of the ionic charge distribution 
can trigger electro-convection |l4|, |l5], |(| ; b) density in- 
homogeneities due to concentration changes at the elec- 
trodes induce lar ge s cale gravity-driven convection rolls 

0,0HH0,|||||g. 

However there is another potential source of density 
changes: due to the small cell volume of typically < 
lcm 3 , the dissipated electric energy can be the source 
of a significant heating. If the cell were completely ther- 
mally insulated, a solution exposed to an electrical power 
of 500 mW would start to boil after 500 s. While ther- 
mal conduction will confine the overall temperature in- 
crease to smaller values, considerable temperature gradi- 
ents might arise and generate density driven convection. 

Thermally induced convection rolls have been thor- 
oughly studied in thin-layer geometries heated from be- 
low (for a recent review see [|l]) and the side |^2|, |2jj, 
p4| , p5| , p6[ . In contrast, the role of thermal effects was 
previously not examined in electrodeposition. In order 
to quantify the possible temperature gradients, a high 
spatial resolution of the temperature field is necessary. 
We present here measurements of the temperature field 
at both electrodes performed by use of an infrared cam- 
era. The so determined evolution of the temperature field 



can be related to the development of the concentration 
field, which is known from interferometric measurements 
0, 27 . The remainder of the paper is organized as fol- 



lows: Section II describes our experimental setup, sec- 
tion III presents the measurement results, which will be 



discussed in sec. IV 
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II. EXPERIMENTAL SETUP 

The measurements are performed with the infrared 
camera Varioscan 3021-ST from InfraTec, which contains 
a Stirling-cooled HgCdTe-detector with 360 x 240 pixel. 
Using its macro we observe an area of 5.0 x 3.4 cm 2 , 
which yields a spatial resolution of 140 fiia. The max- 
imal image capturing frequency is 1.1 Hz, the thermal 
resolution ± 30 mK. 

While the sensor is sensitive for wavelengths in the 
range of 8-12 /im, glass plates, which are normally used 
as top and bottom plate of the cell, are opaque in this 
region. Therefore the upper cell cover was realized using 
a polyethylene foil stretched over an aluminium frame. 
Due to the flexibility of the polyethylene foil and a small 
overpressure necessary to fill the cell, the plate separation 
of about 0.5 mm is not well defined. The bottom plate 
of the cell consists of a block of Teflon, because this ma- 
terial has a low reflectivity in the infrared. This reduces 
the so called narcissism: the response of the cooled detec- 
tor to its own mirror image. The electrodes are parallel 
zinc wires (Goodfellow 99.99%) of 0.25 mm diameter and 
separated by a distance of 4 cm. Fig.|| shows an image of 
the cell during an experiment. 

The cell is filled with an 0.1 M ZnSC-4 solution pre- 
pared from Merck p. a. chemicals in nondeaerated ultra- 
pure H2O. The measurements are performed at a con- 
stant potential of 20 ± 0.003 V. Due to the current in- 
crease during the electrodeposition process, the average 
electrical power feed Qei increases from 470 mW at the 
beginning of the experiment to 650 mW after 500 s. 
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FIG. 1: Cell used for the infrared measurements. The bottom 
(white) consists of Teflon, the upper cover is a polyethylene 
foil stretched over an aluminium frame. The electrodes are 
parallel zinc wires with a distance of 4 cm. The deposit at the 
cathode has grown for 410 s. 



Because of the modified cell construction, the question 
of comparability with experiments performed in standard 
electrodeposition cells could be raised. Therefore we cal- 
culate the heat-flux kiAi for the confining plates, where 
Ai and ki are the area and the heat transfer coefficient 
of the plate, ft is important to keep in mind, that 1/fcj 
is equivalent to l/a% + l/a.2 + Azi/A.;. Here ot\ and a.^ 
are the heat transition numbers from solution to plate 
resp. plate to air, \ is the heat conductivity and Azi the 
thickness of the plate. 

" 2i 



Inserting the material parameters of our setup 1 28 



and using only the area between the two electrodes, we 
derive a kiAi of 16.2 mW per Kelvin temperature dif- 
ference for the polyethylene foil and 11.5 mW/K for the 
Teflon plate, while a typical glass plate (Schott BK 7, 
6.3 mm thick) yields 16mW/K. So in a first approxima- 
tion our setup is thermally equivalent to a standard elec- 
trodeposition cell. 



III. EXPERIMENTAL RESULTS 

Fig. Q shows the growing deposit at the cathode, which 
belongs to the homogenous morphology |3(| . It is charac- 
terized by tip-splitting but retaining a growing front par- 
allel to the cathode. After 500 s a Hecker transition Q 
takes place and the growth front breaks up into more 
spatially localized zones of active development. 

Immediately before each experiment the infrared cam- 
era takes a zero image, which is then subtracted from all 
images taken during the experiment. Therefore the ther- 
mographies solely depict the temperature increase. Fig. || 
gives an example of such a thermography after 280 s, the 
scale at the right describes the temperature increase with 
respect to the beginning of the experiment. The white 



line marks the position of the anodic zinc wire. The 
temperature decrease at the left hand side of the image 
is caused by the thermal conductivity of the aluminium 
frame, the warm "island" in the middle is due to the 
inhomogeneity of the cell thickness. 




FIG. 2: Thermography taken 280 s after the start of the ex- 
periment. The size of the image is 5.0 x 3.4 cm 2 . The white 
line corresponds to the position of the anodic zinc wire, the 
rows between the black lines are averaged to produce Fig. [| 



In order to improve the signal to noise ratio, zones of 
spatial homogeneity and a width of 9.1mm are chosen 
by visual inspection. Such a zone is depicted in Fig. ^ 
with two parallel black lines. Inside this zone all rows 
are averaged, yielding a temperature increase AT(j/,t), 
which is only a function of distance to the cathode y and 
time. 

Fig. H shows the evolution of AT in the neighborhood 
of the anode. It is clearly visible that the temperature 
increase at the anode itself lags behind with respect to 
the one observed in the middle of the cell. 

Fig. [| illustrates the evolution of the temperature field 
at the cathode. The most prominent feature is the ex- 
istence of a local temperature maximum, denoted with 
small arrows. This maximum moves towards the middle 
of the cell, where a plateau of spatially constant temper- 
ature is located. 

In order to characterize the heating process in the cell, 
we pick the temperature at a distance of 25.6 mm to the 
cathode (this is approximately the location, where the 
two convection rolls emerging from the electrodes finally 
meet, as will be discussed in Section [V B| ). Referring to 
this point, we will speak of the bulk in the following. In 
Fig. | this temperature increase ATbuik is given for the 
two experiments presented in Fig. || and [|. The fact that 
the two measurements are almost identical reflects the 
reproducibility of the experiment. The solid line is a fit 
to an exponential function, which is motivated in detail 
in Section IV C. 

The arrows in Fig. || represent the location, where the 
temperature is 20 mK smaller than ATbuik- The distance 
of these points with respect to the anode is denoted L a (t) 
and is shown in Fig. |^ (a) . Its monotonous increase with 
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FIG. 3: Temporal evolution of the temperature field at the 
anode, which is located at 40 mm. All data are averaged over 
a width of 9.1mm. The arrows indicate the position where 
the deviation of AT from the bulk temperature becomes less 
than 20 mK. 



FIG. 5: Temperature increase in the bulk of the cell at a 
distance of 14.4 mm to the anode and 25.6 mm to the cathode. 
The o correspond to the experiment presented in Fig. ^ the 
■ to the one in Fig. [| The solid line is a fit of Equation 5] to 
the o. 
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FIG. 4: Temporal evolution of the temperature field at the 
cathode. All data are averaged over a width of 9.1mm. The 
arrows mark the position of the temperature maximum. 



time is fit ted by a power law, which will be explained 
in Section IV A . Correspondingly, Fig. || (b) shows the 
growth of the distance L c (t) between the the location of 
the temperature maximum in Fig. ^ and the cathode. 
The straight line is a fit t o all data with t > 100 s, which 



will be motivated in Sec. IV B 



IV. DISCUSSION AND CONCLUSIONS 

In principle, the energy balance involves three con- 
tributions: the electrical energy feed into the cell, the 
dissipated ohmic heat and the chemical reaction energy. 
However the short calculation presented in the appendix 
supports the assumption that the chemical energy con- 



tribution is mostly irrelevant. 

The ohmic heat dissipated at some position in the cell 
will be proportional to the local resistivity p(y, t) in a 
one-dimensional model, while p will depend on the local 
ion concentration c(y,t). In the next three subsections 
we will compare the evolution of the temperature field 
at the anode, at the cathode and in the bulk with the 
development of the concentration field, which is known 
from interferometric measurements . 



A. Evolution of the temperature field at the anode 

At the anode Zn 2+ ions go into solution, increasing 
the local concentration and therefore density. While this 
denser solution sinks down to the bottom plate, it gets 
replaced by less dense bulk solution. This mechanism 
drives a convection roll of size L 



16, 19 . According 



to the fluid dynamical description there are two growth 
regimes: Initially during the so called immiscible fluid 
regime L will grow with t , while after some time there 
will be a crossover to the diffusion hindered spreading 
regime with a t 5 growth law. Chazalviel et al. |ll| 
showed, that both the fluid velocity v and the ion con- 
centration decrease with increasing distance to the elec- 
trode. The distance, where v has dropped to zero and 
coincidently c has fallen to Cbuik, defines L. 

Figure |^ shows that the temperature increase in the 
region between the anode and the arrows lags behind the 
one observed in the bulk. As c > Cbuik translates into 
reduced ohmic heating, we identify the position of the 
arrows with the end of the anodic convection roll. A fit 
to its length with: 



L a (t) = at b 



(1) 
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FIG. 6: Temporal evolution of: (a) the distance between the 
anode and the point, where the temperature starts to deviate 
from the bulk. The solid line is a fit with Equation (b) the 
distance between the temperature maximum and the initial 
cathode position. The fit is performed with Equation |^ to the 
data indicated with filled symbols. 



we perform a linear fit with: 

L c (t) = v c t + l c (2) 

which is shown in Fig. ^ (b). We derive v c = 21.2 ± 
0.4 yLtm/s, which agrees well with 20.7 ± 0.8 /zm/s front 
velocity determined from photographs of the deposit. 

l c is found to be 2 ± 0.1 mm. This finite distance can 
be explained by the absence of heat production in the 
metallic deposit because of its low resistivity and the fact 
that its heat conductivity is 190 times higher than water. 
So the deposit is an effective heat sink, the resulting heat 
flux shifts the temperature maximum into the cell. 

If the cathodic convection roll, apart from the fact that 
it starts at the actual front of the deposit, grows in the 
same way as the anodic roll, they meet 450 s after the 
beginning of the experiment at a distance of 25 mm to 
the cathode. This corresponds to the observed change in 
morphology after that time. 



C. Temperature evolution in the bulk 

The temperature increase observed in Fig. [5] can be 
modeled, if we assume, that the whole cell shares the 
constant bulk properties c and p and therefore T. The 
supplied electrical power Q c \ would then be compensated 
by the heating of the system with heat capacity C and 
the heat flow Qrow- 

Qfiaw = -(T-T )J2hAi (3) 

i 

where To represents the ambient temperature and T the 
temperature inside the electrolyte. If we assume Q c \ to 
be constant, the corresponding differential equation: 



is presented in Fig ^ (a). It yields: a = 0.13 ± 0.01 mm 
and b = 0.78 ± 0.01, which indicates that the convection 
roll is in the immiscible fluid regime for the whole run of 
the experiment. According to the scaling analysis pre- 
sented in Ref. p3| , this corresponds to an average plate 
separation of about 650 /xm. 

B. Evolution of the temperature field at the 
cathode 

Due to the growing deposit, there is a zone of ion deple- 
tion in the vicinity of the cathode, the size of which will 
be affected by the convection roll driven by the occur- 
ing density difference. As the decreased c leads to higher 
dissipation, the increased temperature denoted by the ar- 
rows in Fig. |4] is qualitatively explained. The distance of 
the arrows to the initial cathode position L c (t) should 
correspond to the actual size of the deposit. Therefore 



Qe^iT-To^hAi + C— (4) 

2 

has the straightforward solution: 

(T - T ) = T final (l - e -V T ) (5) 

Here Tfj na i = Q c \/ J^. kiAi denotes the finally reached 
temperature difference and t = C / '^2 i kiAi is the time 
constant of the heating up. A fit of liquation || to the ex- 
perimental data is displayed in Fig. g. It yields: Tfi na i = 
9 ± 0.2 K. This translates to an overall heat-flux J2i kiAi 
of 60 mW per K temperature difference. A comparison 
of this result with the values of kiAi calculated in Sec. II 
shows that about 50% of the heat flux takes place through 
the top and the bottom plate of the cell. The heat flux 
through the side walls, the electrodes and the plate area 
beyond the electrodes accounts for the rest. 
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D. Influence of the temperature gradients on the 
convection rolls 

In order to estimate the influence of the temperature 
gradients on the concentration driven convection cur- 
rents, we plot in Fig. [7] the difference between ATbuik 
and a) the cathodic temperature maximum and b) the 
anode. These temperature gradients result in a down- 
flow at the anode and an up-flow at the growth front, so 
they act as an additional driving. 
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FIG. 7: Temperature difference between the bulk and o: the 
cathodic temperature maximum and ■: the temperature at 
the anode. 
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FIG. 8: Temperature dependency of the densities at the anode 
and in the bulk. The • correspond to c = 1.0 m ZnSd, the 
o to Cbuik = 0.1 m ZnSOzt. 



For a more quantitative determination of their con- 
tribution, the temperature dependency of the density 
was measured for different concentrations. The applied 
density measurement instrument DMA 5000 from Anton 
Paar has an accuracy better than 50/ig/cm 3 . 



After 400 s the temperature of the 0.1 M ZnS0 4 bulk 
solution has reached about 30 C which corresponds to 
a density of 1012.3 mg/cm 3 . In accordance with the 
measurements presented in Ref. @ and the theory in 
Rcf. fill , we estimate c at the anode for that time as 
0.25 M. This corresponds to a density of 1036.4 mg/cm 3 
for T = 30° C and 1037.4 mg/cm 3 for the actually mea- 
sured T w 27° C. So the contribution of the temperature 
gradient to the overall density difference at the anode is 
about 4% as visualized in Fig. ||. 

At the cathode c has reached zero at that time, p of 
H 2 is 995.68 mg/cm 3 for T = 30° C and 995.38 mg/cm 3 
for the T — 31° C at the maximum. This results in a 
2% contribution of the temperature gradient to the total 
density contrast. 

As our applied potential is above average for standard 
electrodeposition experiments, we conclude that temper- 
ature inhomogeneities will only weakly contribute to the 
density driven convection rolls. This result justifies with 
hindsight the use of the theoretical description of Chaz- 
alviel et al.|Q. Finally it should be remarked that 
morphologies like stringy pjl ], where the zone of active 
growth is restricted to few small spots with very high 
local current densities, may differ substantially from our 
results. 
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V. APPENDIX 

The standard enthalpy of formation A^H Q of Zn 2+ 
ions in an infinitely diluted solutions is — 153.89kJ/mol 
p2|. It contains three different contributions: 1) the en- 
ergy necessary to liberate an atom from the surrounding 
lattice, 2) the energy needed to ionize the atom and 3) 
the hydration energy, which is set free when the water 
dipoles surround the ion. Only the last term depends 
on the concentration of the solution, decreasing with the 
number of ions already dissolved. 

As the reaction rate at the electrodes is directly pro- 
portional to the electrical current /, so is the chemical 
power Qchcm : 

n A b g9 r m 

Vchcm — — 7^ i (V) 

r z 

where F is the Faraday constant (96485 As/mol) and z 
the charge number. For zinc we derive Q c hem/I = - 
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0.8mW/mA. In the vicinity of the anode this number 
leads to additional heating due to the transformation of 
chemical energy. At the cathode the chemical energy 
stored in the newly produced zinc has to be subtracted 
from the overall heat production, but does not lead to a 
direct cooling of the cathode. 

With our experimental conditions an average current 
of 30 mA feeds an electrical power of 600 mW into the 
cell. Compared to that Qchcm is 24 mW which is about 
4% of the electrical power. 

At the cathode the ion concentration drops fast to zero, 
which justifies the approximation of an infinitely diluted 



solution. At the anode the ion concentration increases 
during the whole run of the experiment. Therefore the 
hydration energy of the newly produced Zn 2+ ions de- 
creases and in consequence the heat, which is released 
from the chemical reaction decrease as well. So the over- 
all heat production from chemical energy at the anode 
is smaller than estimated here and therefore mostly ir- 
relevant for the energy balance. This conclusion agrees 
with the fact that the anode is the coldest point of the 
cell as clearly shown in Fig in spite of this extra heat 
production there. 
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